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FIELP QF THE INVENTION 

This invention relates to improved methods for building geologically plausible 
3D velocity models for use in 3D prestack seismic imaging, in particular for use in 
curved-ray prestack time migration and in such migration followed by prestack depth 
10 migration, and including utilizing geologic data and geostatistical techniques for 
integrating velocity information derived from seismic data, so called "soft data", with 
well data, or so called "hard data", such as checkshot surveys and sonic logs. 



BACKGROUND OF THE INVENTION 

15 Introduction 

In seismic exploration seismic reflection data from subsurface layers are 
collected by multiple receivers typically at the earth's surface. The more recent 3D 
seismic acquisitions cover large areas, collecting extensive 3D seismic data for use in 
subsurface imaging. 

20 Seismic "migration" corrects and improves initial assumptions of near 

horizontal layering in an attempt to model geophysical realities such as dips, 
discontinuities, and curvature of formations. Seismic migration is typically the 
culmination of the image processing for the seismic data, with the goal of producing 
detailed pictures for use in the interpretation of subsurface geologic structures. Such 

25 detailed pictures are important for prospect generation and reservoir characterization, 
such as lithology, fluid prediction, and pore pressure prediction, as well as for reservoir 
volume estimation. 

3D seismic interpretation of migrated data aspires to produce an accurate 
mapping of the subsurface structures, important for oil and gas exploration. Such 
30 seismic interpretation profits from an improved focusing and positioning of subsurface 
reflectors in the migrated data. It is one disclosure of the instant invention that such 
improved focusing and positioning can be made possible by migrating the data with a 
more realistic velocity model. 
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Seismic migration of 3D prestack seismic data is practiced in two forms, time 
migration and depth migration. Each requires a 3D velocity model, a description of 
acoustic velocity structures. Time migration is simpler, less resource intensive and 
considered to be less accurate and less sophisticated. Historically, the industry has not 
5 been bothered to build a sophisticated geologically plausible velocity model prior to 
prestack time migration, relying instead on simple assumptions deemed to be fitting the 
technique regarded as less accurate. For accuracy, the industry currently looks to time 
and resource intensive "depth migration", a technique utilizing massive intentions to 
converge on solutions. 

10 One disclosure of the instant invention is that, surprisingly, an effective and 

efficient focusing and positioning of subsurface reflectors is achievable by seismic time 
migration, but such results are highly dependant on the accuracy of the 3D velocity 
model employed therein. In particular, a geologically realistic, accurate 3D velocity 
model produces surprisingly good results when used with curved ray prestack time 

15 migration. 

The relatively new 3D seismic depth imaging services, techniques that have 
opened a rapidly growing market in recent years, are at present highly interpretive, 
quite expensive and very complex processing jobs. They cany the promise of 
illuminating complicated oil traps under geologic complications, such as sub-salt 

20 prospects in deep water so that to economize on the high cost of drilling in deep water, 
such better seismic imaging technology is in demand. However, reducing the cost and 
time of such complex processing is therefore also important. One disclosure of the 
instant invention is a methodology that so reduces the cost and time. 

Interpreters are historically given the task of tying well markers to 3D seismic 

25 data. Calibrating velocity models derived from seismic data and used for migration by 
tying into well information (hard data) has been historically deemed the interpreter's 
task, preformed post migration. Inconsistencies between processing velocities and the 
depth tying velocities, however, create formidable challenges at the interpretation 
stage. It is a significant improvement of the instant invention to have previously dealt 

30 with this calibration issue, prior to and/or during migration. 

The instant invention provides a better solution, to both imaging and tying well 
information. More focused seismic data, produced as a result of early velocity 
calibration with hard data and geologic data and the use of geostatisticalty sensitive 
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trend fitting (together referred to as "iDEPTHing") prior to migration, as taught 
herein, restores better signal-to-noise ratio for deeper exploration. By utilizing the 
methods of the instant invention curved-ray prestack time migration provides excellent 
migrated seismic data quickly to interpreters* and provides superior preliminary data 

5 for constructing a 3D velocity model for prestack depth migration. 
Summary of State of the Industry 

The cunent industry practice of prestack depth migration is best described by 
John W. C. Sherwood (Sherwood, 1989); "How can we routinely, efficiently, and 
economically extract interval velocities with the accuracy need for depth inversion? ... 

10 First, an educational estimate of the interval velocities is made. This is used to create, 
from the regular stacked (may be migrated) section, an approximate depth-interval 
velocity model. This first guess is going to be wrong. But the data can be processed on 
the basis of this model using a full pre- stack depth migration. Then the data are 
examined and a determination is made of how well it is focused. If they don't appear to 

15 be focused properly, the model can be adjusted, both the depths and the interval 
velocities, to bring the data into better focus. The adjustment is repeated until the 
results are satisfactory. The procedure is today being used very intensively by some 
companies." An abundance of articles on the subject of iterative prestack depth 
migration have been published for the last 20 years (Wang et. al., 1991; Lee and 

20 Zhang, 1992; MacKay and Abma, 1992; Liu, 1997). 

U. S. Pat. No. 4,992,996, entitled "Interval velocity analysis and depth 
migration using common reflection point gathers" relates to a method of performing 
velocity analysis while eliminating the effects on weak signals caused by stronger 
signals with an assumption that assumes that the initial velocity model is reasonably 

25 correct and geologically plausible. The current invention provides a method for 
providing such a needed geologically plausible initial velocity model. 

The iterative technique of depth migration often fails to converge, and it often 
foils to provide prospects with its economics within reasonable ranges of error for 
drilling. The iterative technique is based on the assumption that all iterations will 

30 converge to a unique solution most of the times. However, these iterations will 
converge to a geologically acceptable solution only when the initial model is 
geologically plausible. And, as Sherwood stated, the initial model is going to be wrong 
most of the times: the target depth is often wrong based on the result of iteration. 
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Again, the current invention provides a method for providing such a needed 
geologically plausible initial velocity model. 

Another favorite assumption is that imaging velocities are different from 
velocities for depth conversion due to anisotropy. It has been assumed that seismic 
5 velocities are in general 2-13% faster than the well average velocities (Guzman et. Al., 
1997). However, considering errors in seismic velocity measurements 5-10%, 
anisotropy cannot be so easily assumed. Based on hyperbolic moveout, vertical 
velocity gradients and velocity trend reversal will be a significant source of errors. 
Small errors in stacking velocities will be amplified in converting to interval velocities. 
10 As a reference for geophysical solutions being non-unique based on 

geophysical data collected on earth surface, Al-Chalabi analyzed explicitly the non- 
uniqueness in velocity-depth curve in his article (Al-Chalabi, 1997). He pointed out 
that the errors in measuring velocities from semblance clouds make the problem even 
worse. However, the current industry practice assumes the existence of a unique or 
15 true solution. The current invention deals with velocity calibration in an improved 
manner in order to reduce inaccuracies in seismic velocities and non-uniqueness. 

U. S. Pat. No. 5,089,994, emitted 'Tomographic estimation of seismic 
transmission velocities from constant offset depth migrations" relates to a method for 
improving velocity models so that constant-offset migrations estimate consistent 
20 positions for reflectors, including tomographic estimations of seismic transmission 
velocities from constant-offset depth migrations. This costly method has a fundamental 
problem of non-uniqueness between layer velocities and layer boundary positions. The 
current invention is a prerequisite to reasonable tomographic estimation. 

U. S. Pat. No. 5,513,150, entitled "Method of determining 3-D acoustic 
25 velocities for seismic surveys" relates to a method of producing a velocity volume for a 
seismic survey volume, based on two-way time seismic data and process velocity data. 
This method provides a hand-on interactive tool for velocity editing but does not 
recognize the importance of velocity calibration. 

Constructing improved velocity models requires all different sources of velocity 
30 data. Inaccurate seismic velocity (soft) data are abundant and reliable well data (hard) 
data are sparse. Most of the times, inconsistencies between seismic and well data are 
observed. Mostly the inconsistencies between different sources of velocity data are 
attributed to some unexplained physical mechanism. However, it is difficult to access 



4 



WO 2005/010797 



PCT/US2004/023456 



any physical mechanism due to inaccuracies in seismic velocity measurements. The 
following patent stated that the inconsistencies might come from anisotropy. 

U. S. Pat. No. 6,253,157, entitled "Method for efficient manual inversion of 
seismic velocity information" relates to a method of calculating seismic velocity for 

5 migration purposes as a function of subsurface spatial position that gives the seismic 
processing analyst direct control of the resulting migration velocity model This 
method recognizes the different sources of velocity information and the strengths of 
one source often offset the weaknesses of one another. It states that different velocity 
information sources usually give inconsistent estimates of the instantaneous velocity 

10 and interpreted that this inconsistency may be due to anisotropy that is being ignored. 

Anisotropy is difficult to measure because it can only be measured in the 
laboratories. Laboratory measurements indicating that there were significant intrinsic 
anisotropy was observed in West Africa, and anisotropy migration is required for both 
improved imaging and accurate positioning in such cases. 

15 U. S. Pats. No. 5,696,735 and 6,002,642, entitled "Seismic migration using 

offset checkshot data" relates to a method of migrating seismic data using offset 
checkshot survey measurements. This method uses direct travel time measurement 
from offset checkshot survey instead of a 3D velocity model. Although, such offset 
checkshot surveys are not commonly available, they may be advantageously used in 

20 special situations for shear wave imaging of dipping reflectors. 
Development of Instant Invention 

It was postulated that velocity models constructed according to current 
industry standards could be too smooth for depth imaging and sometimes wipe out key 
prospects. Excessive smoothing in velocity modeling could be a consequence of 

25 embedded velocity errors. Often layer boundaries were interpreted based on over- 
migrated or under-migrated results. It was detrimental to include erroneous boundaries 
in a velocity model. 

A variation of the collocated Co-Kriging method was early on tested by the 
instant inventor on post-stack seismic data for modeling geologic features such as geo* 
30 pressure zones. The collocated Co-Kriging method showed that it could incorporate 
big local anomalies with smooth transition boundaries (Lee and Xu, 2000), integrating 
seismic (soft) data and well (hard) data by the zone of influence according to variances. 
Calibration was used in the test to tie wells fbT post-stack depth migration. Focusing 
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or residual velocity or focusing errors could not be determined with post-stack data. 
The primary purpose of the post-stack depth migration was to correct positioning 
errors. 

In a recent subsequent test, in accordance with the instant invention, velocity 
5 calibration to hard and geologic data and geologically sensitive trend fitting (together 
referred to as "iDEPTHing") was incorporated in constructing a geologically plausible 
velocity model for prestack time/depth migration. The test showed that curved-ray 
prestack time migration based on such velocity model gave astonishingly ("eureka") 
good results in terms of imaging steeply dipping seismic events and tying well markers 
10 (Kenney and Lee, 2003; Lee, 2003). Further, the output of such time migration 
provides a superior, cost effective and efficient input for any subsequent depth 
migration. 
Current Invention 

The current invention discloses improved methods for constructing a 

15 premigration velocity model by the preferred steps of: 1) editing to avoid embedding 
velocity errors and to improve lateral velocity trends of seismic velocities; 2) better 
fitting of vertical velocity trends by calibrating seismic (soft) data to well (hard) data 
(checkshot and/or sonic logs); 3) more accurately interpolating velocities and 
calibration scale factors using geostatistical techniques; 4) incorporating geologic 

20 features, such as known stratigraphic horizons, into a stratigraphic grid, when their 
positioning accuracy is verified, and otherwise utilizing known geologic features for 
editing; and 5) avoiding excessive smoothing by using a variation of geostatistical 
collocated co-Kriging for localized anomolies. 

The invention teaches that curved-ray prestack time migration will correctly 

25 accommodate a vertical velocity gradient. Further, the prestack curved ray time 
migration results using the instant invention provide superior focusing of seismic 
events, including steeply dipping events, and excellent well ties. Due to well ties, more 
accurate stratigraphic horizons can be interpreted. A second velocity calibration (to 
hard data and geologic data) and (geologically sensitive) trend fitting (together referred 

30 to as "©EFTMing*') before subsequent prestack depth migration can significantly 
accommodate lateral velocity variations above a prospect. For sub-salt prospects, the 
inventive methods can provide valuable calibrated velocity data in above top salt for 
constructing a geologically plausible initial model for prestack depth migration. 
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Geostatistical methods have been used in other disciplines, sucn as reservoir 
characterization. Interpreters have used velocity calibration. Some predecessors have 
attempted to calibrate well (hard) and seismic (soft) data prior to migration. The 
current invention, however, is the first to specifically disclose a unique (referred to as 
5 "iDEPTHing") process for calibration prior to prestack time and/or depth migration as 
well as its surprising value when combined with curved-ray prestack time migration, 
and when then subsequently used with depth migration. 

SUMMARY OF THE INVENTION 

The instant invention discloses an improved method for constructing a 3D 

10 geologically plausible velocity model that has particular value for efficient and accurate 
prestack imaging. Embodiments of the invention provide: (1) a method for calibrating 
seismic velocity data, appropriately and effectively taking into account hard data, 
geological features and utilizing geostatistically sensitive trend fitting (together referred 
to as "iDEPTHing") for interpolating, which is particularly valuable for producing 

15 RMS velocities for use with curved-ray prestack time migration; (2) a similar method 
for calibrating and trend fitting ("iDEPTHing") interval velocities before prestack 
depth migration, again appropriately and effectively taking into account well (hard) and 
seismic (soft) data as well as geological features and geostatistical interpolation; and 
(3) a method for efficient sequential use of iDEPTHing in curved-ray prestack time 

20 migration followed by further iDEPTHing and prestack depth migration. Advantages 
of the embodiments include providing a quick turnaround for prestack time and depth 
migration images to interpreters and reducing resource-intensive interpretation efforts 
for 3D seismic data. The invention has significant implications for improving oil and 
gas exploration and production technologies, including pore pressure prediction, 

25 prospect evaluation and seismic attribute analysis. 

Proper calibration and trend fitting of seismic RMS and interval velocities before 
prestack migration is a key aspect of the instant inventive method, leading to efficient 
and successful prestack imaging and interpretation of 3D seismic data. Experiment 
shows that a proper velocity calibration and trend fitting prior to curved-ray prestack 

30 time migration can provide seismic data with excellent imaging and accurate well ties. 
Further velocity calibration and trend fitting of interval velocities can provide a further 
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improved geologically plausible velocity model to entrance pre-stack depth migration, 
both in terms of the speed and the quality of the results. 

Editing seismic velocities using geologic constraints is valuable for effective 
calibration. Use of stratigraphic surfaces can further play a helpful role in more 
accurately calibrating seismic velocities. Geostatistical methods (variogram modeling 
and Kriging) are key factors in interpolating seismic velocities and/or calibration scale 
factors. A scale of geological heterogeneities is better used for smoothing instead of an 
arbitrary scale. 

The inventive procedure provides a method of velocity calibration and trend 
fitting (together referred to herein as "iDEPTHing") for velocity modeling prior to 
curved-ray prestack time migration. Such time migration, in turn, offers an improved 
stepping-stone for efficiently and effectively constructing an accurate 3D velocity 
model for prestack depth migration. 

The invention can be implemented in specific preferred procedures for editing, 
calibrating and trend fitting, thereby effectively integrating well (hard) data, geological 
data and seismic (soft) data. The ("iDEPTHing*') procedure can be supplemented by 
interactive software packages. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 A better understanding of the present invention can be obtained when the 

following detailed description of the preferred embodiments are considered in 
conjunction with the following drawings, in which: 

Figure 1: illustrates two basic steps of one embodiment of the current 
inventive method for velocity calibration and trend fitting for prestack time migration 
25 and geologically plausible velocity modeling for prestack depth migration. 

Figure 2: illustrates preferred steps for velocity calibration and trend fitting 
("iDEPTHing") for curved-ray prestack time migration. 

Figure 3: illustrates preferred steps for velocity calibration and trend fitting 
("iDEPTHing"), and for constructing a geologically plausible velocity model for 
30 prestack depth migration. 

Figure 4: illustrates seismic Velocity Types, (a) Stacking (or NMO) velocities 
(red and blue lines) are interpreted from semblance clouds (ellipses) which are 
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generated from seismic data according to velocity analysis (note: personal bias can 
interpret semblance clouds differently); (b) RMS velocities are inferred from stacking 
velocities assuming the earth flat; (c) average velocities by time averaging of interval 
velocities; (d) interval velocities from the Dix equation as applied to prior RMS 
5 velocities. Three neighboring velocities are shown in Figures (b), (c) and (d). 

Figure 5: illustrates Velocity Editing, (a) A base map and velocity locations 
(red cross) comprising selected locations with pre-existing salt contours (yellow and 
blue) and faults (green and black); (b) interval velocities (red lines) with a geological 
constraint envelope (blue) which is computed and checked with geology; (c) time- 
10 interval velocity data from two selected locations (blue and black lines) before editing; 
(d) the same time-interval velocity data from two selected locations after editing 
including resampling, median and damped-Ieast-square filtering. Note: a couple of 
isolated points were moved into the red envelope by applying filtering. 

Figure 6: illustrates velocity editing, (a) A base map showing lateral velocity 
15 trends in colors (cool colors are slow) in three (northern, mid and southern) blocks at 
four seconds; (b) RMS velocities with a geological constraint envelope (blue) for the 
entire area showing velocity trends; (c) southern blocks with slow velocity trends with 
global envelope (green) and local envelope (blue); (d) and (e) mid and northern blocks 
show fester velocity trends. 
20 Figure 7: illustrates geologic markers. Velocity trends between certain geologic 

markers were discerned from rock types and ages. The estimation of vertical velocity 
trends can be improved based on known rock types and again editing can be applied by 
choosing filter parameters and sampling rates in different geologic time windows 
according to information about rocks. A sonic (green) curve and four neighboring 
25 seismic (red) velocities are shown. Horizontal lines with name (MIOC, OLIG, ...) 
indicate the ages of rocks. 

Figure 8: illustrates Checkshot Data, (a) Time vs. Depth curve from measured 
data; (b) time vs. average velocity derived from checkshot data; (c) time vs. interval 
velocity derived from checkshot data. 
30 Figure 9: (a) A checkshot average velocity (bright blue) was compared with 

neighboring seismic velocities (blue lines); (b) a checkshot RMS velocity (bright blue) 
was compared with seismic average velocities (blue lines); and (c) a checkshot interval 
velocity (bright blue) was computed with seismic interval velocities (blue lines). 
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Checkshot RMS and interval velocities were derived from checkshot data for 
comparison. 

Figure 10: illustrates Stratigraphic Modeling, (a) Three stratigraphic surfaces; 
(b) a cross-section with layer boundaries; (c) an upper stratigraphic unit with a well; 
and (d) a lower stratigraphic unit with a well are shown. Two stratigraphic units have 
stratigraphic grids and the shape of grids are irregular but conforming with two 
stratigraphic horizons. 

Figure 11: illustrates Calibration Process, (a) Pre-calibration seismic RMS 
velocity contours; (b) post-calibration seismic RMS velocity contours; and (c) 
Calibration Table, showing tabulated velocities before and after calibration and scale 
factors. Well derricks indicate well locations, numbered counterclockwise starting from 
far right, and thin lines indicate locations of seismic velocity analysis. 

Figure 12: (a) A velocity trend map at three seconds with known salt contours 
(blue and yellow lines) and faults (purple and black lines); and (b) interval velocity data 
show vertical trends with colors. A horizontal line indicates time equals three seconds. 
Warm colors indicate fast velocities and cool colors indicate slower velocities using a 
rainbow scale. Purple color indicates no data zone. 

Figure 13: A variogram map on a stratigraphic horizon, shows NE direction as 
the direction of maximum spatial continuity. Colors indicate normalized variogram 
values. Two red dots and an ellipse are variogram modeling tools for finding the 
direction of maximum continuity and computing variogram values in major and minor 
axes of ellipse. 

Figure 14: (a) A base shows the location of Line A to B (red line) and known 
salt contours and well locations; (b) Line A to B (base map) of curved-ray prestack 
time 

migration result shows excellent imaging of steep salt face reflection (inside yellow 
ellipse), faults. It also shows excellent ties to well markers: Dn (green), Bh (red), Ym 
(yellow) and salt (pink) inside yellow circles. Vertical scales are time and horizontal 
scales are line and trace locations. 

Figure 15: illustrates a comparison of two prestack time migration, (a) Previous 
1999 prestack time migration result; (b) current 2003 prestack time migration result 
shows significantly improved imaging of salt flank and faults, wherein the accuracy of 
the location of salt is critical for defining oil compartments; (c) a depth slice shows oil 
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compartment defined by the location of salt and faults; (d) well planning by dotted red 
line indicates a successful drilling path close to the salt face compared to a foiled 
drilling path (black line). Detailed sand columns are indicated by varying colors. 

Figure 16: illustrates curved ray prestack migration results, (a) prestack time 

5 migration without velocity calibration and trend fitting; (b) prestack time migration 
with velocity calibration and trend fitting showed improved images. The above two 
sections show significantly different well markers (green bar: DN, pink bar: salt) ties 
with seismic data (yellow circles). Interpreters agreed that (bisection shows excellent 
ties to well markers as a result of using RMS velocities with calibration and fitting. 

10 Horizontal scales indicated Line and Trace locations and the vertical scale indicates 
Time in seconds. 

Figure 17: illustrates common image point gathers, showing consistent 
prediction of times between different offset traces, which indicate small velocity errors. 
Small residual errors indicate that the velocity model was accurate for imaging. 

15 



20 



25 
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DESCRIPTIO N OF THE PREFERRED EMBODIMENTS 

Discussion of Terms 
Time And Depth Migration 

Seismic migration methods are classified as either depth migration or time 

5 migration. Depth migration honors lateral as well as vertical variations, utilizing 
different assumptions about the physics of seismic wave propagation. Time migration 
solutions ignore wave-field distortions created by lateral variations in seismic 
velocities. However simplistic, time migration is widely used because it is less 
computationally expensive. Curved-ray time migration is a recent advance in prestack 

10 time migration, which is sensitive to vertical velocity variations. However, it is only 
with the relatively recent advent of clustered PC processing that curved-ray prestack 
time migration has become generally available. 
Checkshot Data 

Checkshot is a method of determining an average velocity as a function of 
15 depth by lowering a geophone into a hole and recording energy from shots fired from 
surface shot holes. Checkshot is often run in addition to a sonic log to supply a 
reference time at the base of a casing to check an integrated time. Typical intervals 
between receivers are 250 to 500 feet. Some checkshot wells are vertical and some are 
deviated. 
20 Sonic Logs 

A sonic log measures interval travel time (ITT), which is the reciprocal of 
interval (P-wave) velocity in microseconds per foot. It is frequently used for porosity 
determination by time-average equation. 
Geostatistics (Isaaks And Srivastava. 1989) 

25 Geological and geophysical data come in discrete and incomplete forms. A 

common task for a geophysicist is to fill up the space between data points with 
reasonable estimates. The quantitative description of geological trends is analyzed by 
the variogram process and the interpolation method incorporating such variogram 
trend information is called Kriging. Geophysical data also comes with different 

30 measurement scales. Sonic logs are measured in every half foot, and seismic velocities 
are measured to 100 milliseconds, which corresponds to ~ 300 feet. Data of such 
different accuracies can be integrated in interpolation by a process called Co -Kriging. 



12 



WO 2005/010797 



PCT/US2004/02J456 



Variogram (Isaaks And Srivastava. 1989) 

In geostatistics, a variogram is a measurement of spatial variability. In common 
words, it indicates how data points in a defined space become uncorrected as the 
distance between points increases. 
5 Geostatistical Interpolation 

Geostatistical interpolation refers to Kriging and Co-Kriging incorporating 
variogram trend information. 
Velocity Functions 

The term 'Velocity functions" is used to indicate common velocity parameters 
10 used in the industry, such as RMS velocity, interval velocity, and average velocity. 
The term is intended to include such or equivalent (new or old) or newly developed or 
alternate velocity functions. 
Selected Seismic Data 

"Selected seismic data" indicates a use of less than a whole set of recorded 
1 5 traces from a 3D seismic data set. 
Hard Pata Source? 

"Hard data sources" indicate sources of well data such as checkshot data, sonic 
data and log data. 
Established 

20 Use of the word "established" in reference to stratigraphic horizons is intended 

to indicate stratigraphic horizons that are accepted as relatively well known and/or are 
believed likely. There should be significant or strong substantiation for the existence of 
such horizons. 

Alternative Collocated Co-Kriging 

25 A variation of the collocated co-Kriging method, called Alternative Collocated 

Co-Kriging, is designed to incorporate geologic features causing large velocity 
anomalies. The method produces a smooth boundary with a large anomaly using 
estimation variances. 

In general, see Background Information Table. 

30 Discussion of Invention 

The current invention is a method for: 1) efficiently generating imaging and 
well ties for 3D seismic data using prestack time migration; (2) improved velocity 
calibration to hard data and geological data, including geostatistically sensitive trend 
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fitting, ("iDEPTHing") of RMS velocities for use in curved-ray prestack time 
migration; and (3) constructing an improved 3D geologically plausible velocity model 
for subsequent prestack depth migration. 

The current invention is based on the following concepts. 
5 • Seismic velocities indicate lateral velocity variation faithfully but tend to 

be inaccurate and unreliable for vertical velocities, due to the hyperbolic moveout 
assumption. 

• Discrepancies between seismic velocities and well velocities are not 
necessarily due to anisotropy, as assumed by many people. They may be due to, and 

10 explained by, inaccuracies in seismic velocity analysis. Early calibration and trend 
fitting ('iDEPTHing") can surprisingly resolve many such inconsistencies in velocities 
derived from well (hard) data and seismic (soft) data (Figure 9). The prominent 
discrepancies between seismic (soft) and checkshot (hard) data are due to geopressure 
in Figure 9. The procedure can improve early estimated velocities in RMS, average and 

1 5 interval velocities. 

• An early accurate velocity model requires more than just seismic 

velocities from seismic data. Various other usable sources of velocity data are 

frequently available, and should be appropriately integrated into the model, such as 

checkshot, sonic, log, well markers, proximity surveys, VSP data and other geologic 

20 markers or descriptions. (The latter are typically referred to as hard data.) The value 

# 

of checkshot and tog data in adjacent blocks, in particular, is highlighted by the current 
invention. 

• Surprisingly, curved-ray prestack time migration can generate excellent 
seismic imaging and well ties if one first calibrates and trend fits (preferably edited) 

25 RMS velocities (with preferably edited) available hard data and geophysical data. A 
further improved 3D geologically plausible velocity model subsequently can be 
constructed by interpreting prestack time migration results. 

• Seismic velocities are preferably edited with the use of geologic 
constraints, to the extent available, and in particular, seismic interval velocities should 

30 be adjusted to be stable and be in agreement in terms of predictable rock properties. 
RMS velocities and average velocities are best derived from edited interval velocities. 
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• Automated editing tools are particularly helpful for efficient volume 

editing. 

• Sophisticated calibration tools are particularly useful so that sparse well 
data can be smartly interpolated. 

5 • Geostatistical Interpolation (Kriging and Co-Kriging) is best utilized for 

smart interpolation, honoring different correlation lengths in lateral directions. Kriging 
is valuable for honoring the anisotropic correlation lengths. 

• Geostatistical methods are best used to integrate the well (hard) data 
and seismic (soft) data. An alternative method could use a geological heterogeneity 

10 scale. These methods best replace excessive smoothing, now popular in the industry. 

• A set of velocity tools best enable the current invention to proceed in 
steps, including editing, building stratigraphic units, and geostatistical interpolation 
methods and calibration including trend fitting* For quality control, the calibration and 
trend fitting results can be checked along well tracks and on curved well planes 

15 between wells. 

Discussion of the Figures 

Figure 1 is a schematic flow chart illustrating a preferred embodiment of the 
current invention. In Steps 10-12 indicate velocity calibration of seismic RMS 
velocities with hard data and geostatistically sensitive trend fitting (together sometimes 

20 referred to as "iDEPTHing") before curved-ray prestack time migration, one key 
aspect of the current invention. Curved-ray prestack time migration accommodates 
vertical velocity gradients. The result surprisingly images well all seismic events and 
ties well markers. Such result alone can save interpreters significant time by not having 
to task for well ties. 

25 Subsequently, migrated time data can be converted to depth using an average 

velocity cube, converted from the RMS velocity cube. The depth of any stratigraphic 
horizons produced should be accurate for velocity depth modeling in the next steps. 
Steps 20-22 indicate further velocity calibration, again with hard data and 
geostatistically sensitive trend fitting ("DEPTHing"), for subsequent prestack depth 

30 migration. Such comprises an efficient method for constructing a geologically 
plausible initial interval velocity model. This model calibrated and trend fitted, takes 
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advantage of the further information derived from the migrated time data. New 
stratigraphic horizons may have been identified. 

According to the current industry practice, iterative pre stack depth migration 
relies on cycling through Steps 24-28 without performing Steps 10-14 and 20-22. 

5 According to current industry practice, a macro model is constructed using interpreted 
horizons and interpolating velocity fields. This type of model has the following 
characteristics: 1) velocity errors are embedded due to lack of proper editing; 2) 
velocities are smoothed to hide the embedded errors; 3) the positioning of layer 
boundaries are not verified; and 4) vertical velocity trends are deviated whenever 

10 velocity reversals or large velocity gradients exist. Steps 10-14 and 20-22 can improve 
the initial velocity model by avoiding embedded velocity errors, avoiding excessive 
smoothing and by trend fitting using geostatistical Kriging. 

For preferred embodiments of the present invention, more detailed steps are 
described in subsequent flow charts, indicating in more detail velocity editing, 

15 calibration and trend fitting in Figure 2 and constructing a geologically plausible 
velocity model in Figure 3. 

Figure 2 illustrates one preferred method for a preferred embodiment, 
indicating velocity editing, calibration and trend fitting for use in curved-ray prestack 
time migration. (Step 40) Velocities are compiled from selected 3D prestack seismic 

20 data and interval and RMS velocities are edited honoring geologic trends that exist in 
the survey area. Seismic velocities can be edited using envelopes of vertical trends, 
based on rock properties and/or using computed envelopes. (See Figures 4, 5, and 6). 
Regional envelopes may be needed for an extensive area (Figure 6). Envelopes form 
constraints. Editing can include resampling, to remove anomalies and/or applying 

25 median and damped-least-squares filters to RMS velocities. The RMS velocities and 
interval velocities may be edited with the help of interactive windows. 

Lateral trends can be computed by geostatistical variogram modeling. (See 
Figures 12, 13). Geological maps, faults, geo-pressure zones, geological markers and 
salt intrusion can be viewed on an interactive workstation, including on dual windows, 

30 with one window for interval velocities and one window for RMS velocities. 

More particularly, the following steps can be taken for editing. Referencing 
step 41 of Figure 2, one can convert computed seismic RMS velocities to interval 
velocities using the Dix equation, and to average velocities. Interval velocities are 
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defined at the center of the layers (See Figure 4). One can compute envelope interval 
(unctions by specifying upper and lower limits based on geologic constraints, to the 
extent available. One can then delete or correct erratic functions of picks, lying 
outside of envelope functions. (See Figure 5). Referencing step 42 of Figure 2, an 

5 alternative possibility is to resample the seismic data and/or also to apply median and 
damped-least-square filters on the RMS velocity domain. If rock properties are known 
between certain geologic markers (See Figure 7), one can apply filters for user-defined 
time windows. Time-varying re-sampling is an option to consider for velocity trends 
between certain geologic markers. Any deleted velocity functions could alternately be 

10 assigned lower weights. Anomolous interval velocities can be corrected to their 
envelope value. 

(Step 44) RMS velocities can now be calibrated, preferably using selected 
stratigraphic horizons which control sedimentation. (See Figure 10). The following 
step can be taken for stratigraphic horizons. 

15 Referencing step 43 of Figure 2, known horizons controlling sedimentation 

styles due to compaction, if any, such as the water bottom horizon in the Gulf of 
Mexico, are selected. Such stratigraphic horizons provide conforming surfaces for 
velocity calibration. A conforming stratigraphic unit is built using, for example, two 
stratigraphic surfaces (See Figure 10). 

20 Referencing step 45 of Figure 2, in order to calibrate RMS velocities, 

preferably using key controlling stratigraphic horizons, variogram models for velocity 
trends from seismic velocities are computed. (See Figures 11, 12 and 13). Checkshot 
and/or other hard well data (step 46) is edited to the extent feasible. (See Figure 7). 
RMS and/or interval velocities are computed at well locations from hard data (step 

25 47). (See Figure 8). The next step is to interpolate seismic velocities to the locations 
of hard well data (step 48). One is now in a position to compute scale factors, 
preferably by dividing checkshot RMS velocities at well locations by seismic RMS 
velocities interpolated to the well locations (step 49). Geostatistical Kriging, using the 
variogram model of velocity trends, is used to interpolate the scale factors to the 

30 selected locations of seismic velocities (step SO). Calibrated seismic velocities are 
computed by applying (such as by multiplying) the interpolated scale factors to the 
seismic velocities (step 51). See Figure 11. 
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(Note: Deviated checkshot and checkshot from adjacent blocks, if available, 
can be included for calibration.) 

Again, checkshot and/or other hard data is also preferably first edited, such as 
by reviewing time vs. depth, time vs. average-velocity and time vs. interval-velocity 
5 relationship (See Figure 8). Some checkshot data may be determined not to be used 
because it measures a local anomaly. 

(Step 52) Kriging or Co-Kriging is applied to the calibrated seismic velocities 
at the select locations to geostatistically interpolate the calibrated velocities to all trace 
locations. The process includes fitting calibrated velocity data to computed velocity 
10 trends, consistent with geologic trends, for geostatistically sensitive trend fitting of the 
calibrated data. 

Referencing step S3 of Figure 2, one can further compute a time-slice of 
seismic velocities and check or match lateral velocity trends with geologic maps, if 
available (See Figure 12). 
15 Referencing step 54 of Figure 2, one can further compute a variogram model of 

seismic velocities and compare or match it with "geologic trends" if available (See 
Figure 13). 

Referencing step 55 of Figure 2, calibrated seismic RMS velocities are 
interpolated by the Kriging process on a stratigraphic grid to trace locations for use in 
20 prestack time migration. Calibrated seismic average velocities can also be 
geostatistically interpolated for time-to-depth conversion and depth interpretation. 

Step 60 of Figure 3 illustrates constructing a geologically plausible velocity 
model for subsequent prestack depth migration by the following steps. 

In Step 71, seismic interval velocities are calibrated analogously to the above. 
25 The scale factor between true depth and migrated depth can be added into the scale 
factors computed from checkshot data. Salt entry points and proximity surveys can 
also be used for calibration in a similar manner to other hard data. 

In Step 72, sonic (interval) velocities are resampled and calibrated. Sonic 
velocities can be calibrated with checkshot data to be consistent with calibrated seismic 
30 velocities. 

For trend fitting as well as for use in calibration (steps 81 and 82), prestack 
time migration results can be interpreted and verified against well marker ties. Key 
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stratigraphic horizons can be updated with time-to-depth conversion of time horizons 
(Figure 14). Variogram models can be computed for each stratigraphic unit. 

In Step 83, calibrated seismic interval velocities are interpolated to all trace 
locations on new-updated stratigraphic horizons for prestack depth migration. 
5 In Step 80, similar geo statistical methods can be used to integrate sonic logs to 

seismic velocities, preserving the local effect of the data of the sonic log. 

In Step 81, calibrated sonic (hard) data and seismic interval velocities (soft 
data) can be integrated by a geo statistical method. Alternative Collocated Co-Kriging 
is a preferred collocated co-Kriging method for such integration. Geologic 
10 heterogeneity scales are preferably used instead of an arbitrary smoothing scale (See 
Figures 14 and IS). Such scales are indicated by known anomalous geologic features. 
Further Discussion 

Velocity calibration has traditionally been applied after migration, during 
interpretation. The instant invention teaches geologically sensitive velocity calibration 
is to be more profitably applied before prestack migration. Further, velocity editing is 
important and helpful before velocity calibration. Prestack migration with 
"iDEPTHing" improves both focusing and well ties. 

Velocity trends can be carefully examined particularly in an extensive (ex. 200) 
block area. Several envelope functions can be computed for editing and preserving 
20 velocity trends toward the directions of rapid variation. 

Stratigraphic horizons should be chosen if the boundary controls sedimentary 
layering. A good example is the water bottom in the Gulf of Mexico. Not every 
velocity layer boundary can be a stratigraphic horizon. 

Checkshot data should be carefully verified and checked for any errors 
25 common in checkshot interval velocity. Checkshot calibration fits not only vertical 
trends but also lateral trends using a variogram model. Some checkshot data should 
not be used for calibration, however, if they represent local geology such as a local 
Carbonates seam. Sonic data can be summed for every 100 feet and block shifted to be 
used for calibration. 

30 Geological features such as geo-pressure zones can be handled by the localized 

integration of sonic (hard) data. The integration is local within the geological 
heterogeneity scale using estimation variances. 
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The current invention teaches that a geologically plausible model can be 
constructed based on a starting isotropic earth assumption. Curved-ray prestack time 
migration results will tie well markers and provide a helpful guide for drilling 
prospects. Residual processing will enhance focusing and signal-to-noise ratio (See 
5 Figure 17). 

If the residual moveout becomes significant, anisotropy may be considered. 
Anisotropy should be measured from a laboratory, however, and anisotropic prestack 
migration is needed for further improvement of seismic imaging and well ties. An 
anisotropic factor model can be constructed from residual velocity errors with 
10 calibration with anisotropy data from laboratory measurements. 

The current invention can be developed to utilize shear-wave seismic data or 
multi-component seismic data, with minor modifications. 
HfrtQrigaJ Perspective 

15 Curved-ray prestack time migration has been known in the industry for many 

years but was not practically available to most in the industry until PC clusters became 
a viable hardware solution. In the oil industry as major vendors started moving toward 
PC clusters, new algorithms, such as curved-ray prestack time migration, offering 
further accuracy became, generally, available. 

20 Prestack time migration has typically been used for AVO studies and for 

seismic attribute studies. Traditionally, the industry has not paid attention to editing 
and improving RMS velocities for such time migration because time migration has 
been considered primarily as a process for increasing signal-to-noise ratios by 
collapsing diffractions. Because RMS velocities typically are smooth and 

25 monotonically increasing, and time migration results do not show any significant 
seismic image distortions therein due to velocity errors, it has historically been 
considered acceptable to use RMS velocities in time migration without checking the 
corresponding interval velocities converted from RMS velocities. 

The present inventor is the first to disclose and teach that velocity calibration 

30 and trend fitting for prestack time migration greatly enhances the value of the results, 
for imaging and tying well markers, and can expedite more complicated depth 
migration. In December, 2002, the instant inventor persuaded Mr. Michael Kenney, a 
managing partner of Summit Energy Co. to try an experiment. The result was a 
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remarkable success, correctly imaging all highly dipping salt faces and tying well 
markers on all sides. The results showed everything that Mr. Michael Kenney dreamed 
of seeing, including many interpretable geologic features. 

The experiment substantiated that velocity calibration and trend fitting prior to 
5 curved-ray prestack time migration could yield amazing results in terms of imaging and 
well tying. Further, such prestack time migration could advantageously be used as an 
initial step for prestack depth migration and secure significant economies. Further 
calibration and trend fitting can constrain a velocity-depth model after prestack time 
migration. 

10 The 2QQ2/2QQ? Experiment ("Eureka") 

The experiment tested the value of constructing a 3D geologically plausible 
velocity model prior to curved ray prestack time migration. The reprocessing of a 
proprietary 3D data set attempted to achieve the goal of extracting a high quality salt 

15 flank image as well as tying into wells around the salt. Seismic velocities were edited 
in the interval velocity domain using geologic constraints. RMS velocities were 
calibrated with geostatistically interpolated scale factors based on checkshot data and 
further trend fitted to trace locations. The water bottom was used as a stratigraphic 
horizon. The whole 3D seismic volume was migrated by curved-ray prestack time 

20 migration. The migration results provided excellent well ties on all sides of the salt 
structure in the survey area (Figures 14 and 16). The common image point gathers 
show consistencies in predicting times due to velocity accuracy (Figure 17). Several 
reliable stratigraphic horizons can be interpreted from the curved-ray prestack time 
migration results. The calibration and trend fitting ("iDEPTHing") of interval velocities 

25 included well marker miss-ties and known salt entry points. 

The foregoing description of preferred embodiments of the invention is 
presented for purposes of illustration and description, and is not intended to be 
exhaustive or to limit the invention to the precise form or embodiment disclosed. The 
description was selected to best explain the principles of the invention and their 

30 practical application to enable others skilled in the art to best utilize the invention in 
various embodiments. Various modifications as are best suited to the particular use 
are contemplated. It is intended that the scope of the invention is not to be limited by 
the specification, but to be defined by the claims set forth below. 
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